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Abstract 



The Perona-Malik equation is a celebrated example of forward-backward parabolic 
equation. The forward behavior takes place in the so-called subcritical region, in which 
the gradient of the solution is smaller than a fixed threshold. In this paper we show 
that this subcritical region evolves in a different way in the following three cases: dimen- 
sion one, radial solutions in dimension greater than one, general solutions in dimension 
greater than one. 

In the first case subcritical regions increase, but there is no estimate on the expansion 
rate. In the second case they expand with a positive rate and always spread over the 
whole domain after a finite time, depending only on the (outer) radius of the domain. 
As a by-product, we obtain a non-existence result for global-in-time classical radial 
solutions with large enough gradient. In the third case we show an example where 
subcritical regions do not expand. 

Our proofs exploit comparison principles for suitable degenerate and non-smooth 
free boundary problems. 
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parabolic equation, moving domains, subsolutions and supersolutions. 



1 Introduction 



In this paper we consider the Perona- Malik equation 



u t (x, t) — div 



( 



l + |V«(x,t)| 2 



Vu(x,t) 



) 



V(i,t)efix[0,T), 



(1.1) 



where f2 C M. n is an open set and T > 0. This equation is the formal gradient flow of 
the functional 



The convex-concave behavior of the integrand makes (1.1) a forward-backward par- 
tial differential equation of parabolic type. The forward (or subcritical) region is the set 
of points (x, t) where |Vw(x,t)| < 1, the backward (or supercritical) region is the set of 
points where \Vu(x,t)\ > 1. 

This equation was introduced by P. Perona and J. Malik [17] in the context of 
image denoising. They considered equation (1.1) in a rectangle Q C M 2 , with an initial 
condition u(x, 0) = u (x) representing the grey level of a (noisy) picture, and Neumann 
boundary conditions. For increasing values of t the functions x — > u(x,t) are thought 
as successively restored versions of uo(x). 

The heuristic idea is that small disturbances, corresponding to small values of the 
gradient, are smoothed out because of the diffusion which takes place in the forward 
regions. On the other hand, sharp edges correspond to large values of the gradient and 
therefore they are expected to be enhanced by the backward character of the equation 
in supercritical regions. This phenomenology has been actually observed in numerical 
experiments, which also reveal an unexpected stability (see [7, 8, 14]). This discrepancy 
between the practical efficacy of (1.1) and its analytical ill-posedness has been called 
"Perona-Malik paradox" after S. Kichenassamy's paper [16]. 

In the last fifteen years the paradox has been investigated in numerous papers. Sev- 
eral authors proved well posedness results for approximations of (1.1) obtained via space 
discretization [4, 9] or convolution [5], time delay [1], fractional derivatives [13], fourth 
order regularization [2], simplified nonlinearities [3]. The behavior of such approxima- 
tions as the suitable parameter goes to is a much more challenging problem. As far as 
we know, results in this direction have been obtained only for the semidiscrete scheme 
in dimension one. In this case the authors [9] (see also [4]) proved that approximated 
solutions converge to a limit, and under reasonable assumptions this limit is a classical 
solution of (1.1) inside its subcritical region. 

Despite of these partial progresses, a solution of the paradox is still far away. We 
remind that a solution of the paradox is a notion of weak solution for (1.1) which exists 
for large classes of initial data (for example in BV(Vt) or SBV(£l)), is reasonably stable, 
and to which reasonable approximations converge. 
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In this direction, K. Zhang [19] (see also [6, 18]) showed that the class of Lipschitz 
solutions is far from being a solution to the paradox. Indeed he proved that for any 
nonconstant smooth initial condition, even if subcritical, the Neumann boundary value 
problem admits infinitely many (pathological) Lipschitz solutions. 

Classical solutions (namely solutions which are at least of class C l ) have also been 
investigated in the last decade. B. Kawohl and N. Kutev [15] observed that global- 
in-time classical solutions exist if the initial condition is subcritical, while in [16] it is 
remarked that local-in-time classical solutions cannot exist unless the initial condition 
is very regular in its supercritical region. Moreover the authors proved in [10] that in 
dimension one there exists a dense set of initial data for which the Cauchy problem 
with Neumann boundary conditions admits a local-in-time classical solution of class 
C 2 ' 1 (namely with two continuous derivatives with respect to space variables, and one 
continuous derivative with respect to time). On the other hand, such solutions cannot 
be global if the initial condition has a nonempty supercritical region (see [15] and [12]). 

Quite surprisingly, things are not so drastic in dimension greater than one. Indeed 
the authors proved in [11] that global-in-time radial solutions of class C 2 ' 1 do exist for 
some classes of initial data with nonempty supercritical region. 

For classical solutions one can define the family of open sets 

I-(t) :={x Ett: \Vu(x,t)\ < 1} Vt G [0,T). 

This is the family of subcritical regions. Its behavior as t varies is the object of 
this paper. We point out that this definition is purely local, in the sense that it does 
not depend on the boundary conditions. We show three situations in which subcritical 
regions evolve in a different way. As far as we know, these different behaviors had 
not been explicitly reported in numerical experiments. The reason is maybe that in an 
unstable framework it is always difficult to distinguish what depends on the model itself, 
and what depends on the implementation. In any case, we leave to numerical analysts 
and applied mathematicians any comment about the significance of these results for the 
model and its practical applications. 

The one dimensional case In dimension one subcritical regions don't shrink, namely 

r{s)cr(t) whenever < s < t < T. (1.2) 

This fact had already been proved in [15] under some structural assumptions on the 
initial condition, afterwards removed in [12]. 

In this paper we present an alternative proof (see Theorem 2.1) based on a compari- 
son principle, which only requires C 1 regularity, and works substantially for all equations 
which are the gradient flow of a nonconvex functional. This proof gives us the oppor- 
tunity to show, in a simpler setting, the method which is fundamental in the next case, 
when the result was not known before. 
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This result is optimal in the sense that it may happen that the subcritical region 
is the same for every t G [0, T). An example is provided by the local-in-time solutions 
constructed in [10]. 

The radial case Let us consider a radial solution of (1.1) defined in a ball or an anulus. 
Then (1.2) holds true also in this case. If moreover /~(0) 7^ 0, then the inclusion is 
strict whenever s < t, and there exists T > such that I~(t) = f2 for every t > T . The 
value of T depends on the outer radius of fl, but is independent on u. In other words, 
supercritical and critical regions disappear after a finite time depending only on Q. 

If the solution survives and remains regular up to T , then it becomes subcritical 
and there are no more obstructions to global existence. This is actually what happens 
in the classes of global-in-time radial solutions constructed in [11]. 

In Theorem 2.3 we show that this is not always the case. If the initial condition 
has a supercritical region where the gradient is large enough, and this supercritical 
region is surrounded by subcritical regions, then there is no C l classical solution with 
T > To, independently on the boundary conditions. The reason is that the maximum 
of \Vu(x, t)\ in this supercritical region is a function of time which cannot decrease too 
fast. On the other hand, the supercritical region must disappear after T , and thus this 
maximum has not enough time to decrease from its large initial value up to 1. 

The nonradial case Roughly speaking, in the radial case each interface between the 
subcritical and the supercritical region is a circle which evolves with velocity proportional 
to the inverse of its radius. This reminded us of the mean curvature motion, and in a 
first moment led us to suspect that this interface could evolve in a similar way also for 
nonradial solutions. In particular it seemed reasonable that this interface could evolve 
in such a way that supercritical regions tend to shrink, at least where they are convex. 

In Theorem 2.4 we show that this is not the case. We prove indeed that there exists 
an initial condition u in two variables with the following properties. The supercritical 
region of Mo is convex in a neighborhood of the origin, and any local-in-time solution 
starting from uq, independently on the boundary conditions, has a supercritical region 
which expands in a neighborhood of the origin. In particular any such solution does not 
satisfy (1.2). 

Evolution of supercritical regions If subcritical regions expand, then supercritical re- 
gions shrink. This is true, but not obvious. Indeed it could happen that both the 
subcritical and the supercritical region expand at the expense of the critical region 
where \Wu(x, t)\ = 1. Nevertheless one can prove that this is not the case. It is enough 
to apply the techniques of this paper to the equation obtained by reversing the time. In 
this way the role of subcritical and supercritical regions is just exchanged. 
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Connection with free boundary problems The evolution of subcritical regions is itself 
a free boundary problem. Unfortunately it involves a forward-backward equation. Af- 
ter some variable changes we reduce ourselves to more standard situations. Roughly 
speaking, in the new variable v we end up with equations such as 

Vt = \[v v xx , v t = y/v {v xx + A + lower order terms} , 

where A > 0. These equations are satisfied where v > 0, and we are interested in the 
evolution of the region where v is positive. 

The good news are that these equations are forward parabolic. The bad news are that 
they are degenerate and they involve a nonlinear term which is not Lipschitz continuous. 
This complicates things when using comparison principles. 

For the first equation we show (see Theorem 2.5) that the region where v is positive 
does not shrink with time. This is enough to threat the Perona-Malik equation in 
dimension one. For the second equation we show (see Theorem 2.6) that the region 
where v is positive expands with a positive rate depending on A. This implies our 
conclusions for the radial Perona-Malik equation. 

This paper is organized as follows. In section 2 we state our results for the Perona- 
Malik equation and the related free boundary problems. In section 3 we prove these 
results. 



2 Statements 

Throughout this paper we assume that tp G C°°(R) is an even function, hence in partic- 
ular (fi'(0) = 0. We also assume that 

<p"(a)>0 Vcre[0,l), (2.1) 

= 0, (2.2) 

(p"(a) < V<7 > 1. (2.3) 

These assumptions imply that > and < 0. In some statement we also 

need the stronger assumption 

<p"'(l) < 0. (2.4) 

These assumptions are consistent with the concrete case <p(cr) = 2 _1 log(l + a 2 ) of 
the Perona-Malik equation. We refer to Figure 1 for the typical behavior of <//(<r). 
We consider the following equation 

u t = dwL'(\Vu\)^-Y (2.5) 
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which is the natural generalization of (1.1). We believe and we hope that this generality 
simplifies the presentation, and shows more clearly which properties of the nonlinearity 
are essential in each step. For the sake of generality one could also weaken the regularity 
assumptions on (p (we never consider more than 3 derivatives), replace the threshold 
a — 1 in (2.2) with any a±, and weaken (2.1) and (2.3) by asking only that <//' is positive 
in a left-hand neighborhood of U\ and negative a right-hand neighborhood of U\. 

2.1 Main results 

Let us state our results on the evolution of subcritical regions. The first result concerns 
the one dimensional case, where (2.5) reduces to 

u t = (<p'(u x )) x (2.6) 

This form of the equation is suitable for C 1 solutions, because it involves only first 
order derivatives. When the solution is of class C 2,1 , equation (2.6) can of course be 
rewritten as Ut = <p"(u x )u xx . 

Theorem 2.1 Let if G C°°(R) be a function satisfying (2.1) through (2.3). Let x\ < 
X3 < X4 < Xi and T > be real numbers. 

Let u G C 1 ((xi,x 2 ) x [0,T)) be a function satisfying (2.6) in (xi,x 2 ) x [0, T), and 

Mx,0)|<l \/xe(x 3 ,x 4 ). (2.7) 

Then \u x (x,t)\ < 1 for every (x,t) G (£3, £4) x [0, T). 

We point out that in Theorem 2.1 we don't need assumption (2.4), which was used 
in all previous results of the same type. 

Our second result concerns the radial case. Let r := |x| be the radial variable, and 
let u(r,t) be a radial solution. In this case (1.1) becomes 

Ut = ( v \u r )) r + (n-l)^-, (2.8) 

where n is the space dimension. From now on we assume for simplicity that n = 2. The 
general case is completely analogous. 

Theorem 2.2 Let ip G C°°(R) be a function satisfying (2.1) through (2.4). Let < 
ri < r 3 < r 4 < r 2 and T > be real numbers. Let u G C 1 ((ri, r 2 ) x [0, T)) be a function 
satisfying (2.8) (with n = 2) in (ri,r 2 ) x [0,T) ; and 

K(r,0)|<l \/re(r 3 ,r A ). 

Let k := r^ 1 ^2ip' (l)\ip"' \ , and let us set 

V := {(r, t) G (n, r 2 ) x [0, T) : r 3 - k t < r < r 4 + k t}. (2.9) 

Then \u r (r,t)\ < 1 for every (r,t) G V. 
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In other words, this result says that in the radial case the subcritical region expands 
with a rate which is bounded from below by a positive constant k . Figure 2 shows the 
shape of the set T>. The slope of the slanted lines depends on ho- It is clear that when 
t > ij2 — ri)/fco every nonempty initial subcritical region (r 3 ,r 4 ) has invaded the whole 
interval (r\, r 2 ). 

The third result concerns the nonexistence of global-in-time C 1 radial solutions if the 
gradient of the initial condition is too large. We point out that this result is independent 
on the boundary conditions. 

Theorem 2.3 Let tp G C°°(R) be a function satisfying (2.1) through (2.4). Let < 
fi < r 3 < r 4 < r 5 < r 2 and T > be real numbers, and let u G C 1 ((ri,r 2 ) x [0, T)) be 
a solution of (2.8) (with n = 2) such that 

k(r 3 ,0)|<l, k(r 5 ,0)|<l, u r ( U , 0) > 1 + r2 ^ 2 ~ n) ( 2 - 10 ) 
Then we have that 

r 2 (r 2 -ri) 



T < 



Our last result is a counterexample to the expansion of subcritical regions for non- 
radial solutions. Note that the condition we impose on the initial condition uq depends 
only on the Taylor expansion of Uq of order 3 in a neighborhood of the origin. 

Theorem 2.4 Let y? G C°°(R) be a function such that ip'(l) > and (f"(i) = 0. Let 
Q C ]R 2 be any open set such that (0, 0) G Q. Let u : Q — > R &e any function of class 
C 3 snc/i £/ia£ 

/2 /2 

uo(a;, y) = ^-x + ^-y + hx 2 + k 2 y 2 + h lX 3 + h 2 y 3 + o ((x 2 + y 2 f 2 ) (2.11) 

as (x,y) — > (0,0). Tazs clearly implies that |V-a (0,0)| = 1. Moreover one can choose 
the parameters k±, k 2 , hi, h 2 in such a way that the following properties hold true. 

(1) There exist 5 > 0, a > 0, and a convex function g : (—8, 5) — > (—a, a) such that 

{(x,y) G (-5,5) x (-a, a) : \Vu (x,y)\ > 1} = 

= {(x, y) G (-5, 5) x (-a, a) : y > g(x)} . 

(2) Let T > 0, and let u G C 2 (Q x [0,T)) be a function satisfying (2.5), and the initial 
condition u(x,y,0) = u (x,y) for every (x,y) G f2. Then there exists t G (0, T) 
such that 

|Vu(0,0,t)| > 1 Vte(0,t o ). 

In other words, at time t — the origin lies in the interface which separates the 
subcritical and the (locally convex) supercritical region, while for t G (0,to) the origin 
has been absorbed by the supercritical region. 
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2.2 Heuristics 



In this section we present simple "proofs" of Theorem 2.1 and Theorem 2.2. 

Let us start with Theorem 2.1. Let us assume that u is smooth enough, and that 
one component of the interface between the subcritical and the supercritical region 
is represented by a smooth curve (a(t),t). Just to fix ideas, let us assume that the 
subcritical region lies on the left of the interface, namely where r < a(t). Taking the 
time derivative of the relation u x (a(t),t) = 1 we obtain that 

= u xx (a(t),t)a'(t) +u xt (a(t),t) 

= u xx a'(t) + (p"(u x )u xxx + (p"'(u x )u 2 xx , 

where all the partial derivatives of u in the last line are computed in the point (a(t),t). 
Recalling that (p"(u x (a(t),t)) = 0, we have therefore that 

a' it) = -ip"'(l)u xx (a(t),t). 

Now we have that <p"'(l) < 0, and u xx (a(t),t) > because the subcritical region lies 
on the left of r = a(t). We have thus proved that a'(t) > 0, hence the subcritical region 
tends to expand. 

This "proof" is not rigorous for several reasons: we assumed that u is of class C 3 , 
we assumed that the interface is a smooth curve, we divided by u xx which could be 0. 
Nevertheless we believe that this simple argument is quite explicative. 

Let us consider now Theorem 2.2. As before, we assume that the interface is given 
by a smooth curve (a(t),t), the subcritical region being on the left. Taking the time 
derivative of the relation u r (a(t),t) = 1 we obtain that 

= u rr (a(t), t)a'(t) + u rt (a(t), t) 

a \ in \ mi \ 2 <p"(u r ) tp'(u r ) 

= u rr a (t) + ip {u r )u rrr + if [u r )u rr H j—r u rr r— ^, 

a[t) a z [t) 

where all the partial derivatives of u in the last line are computed in the point (a (£),£). 
Recalling that ip"(u r (a(t),t)) = 0, we have therefore that 

Applying the inequality between the arithmetic and the geometric mean, we thus 
obtain that 



> 2^(1)1^(1)1 > 2^(1)1^(1; 



a(t) r 2 

This "proves" that the subcritical region expands with a rate which is bounded from 
below by a positive constant. The value of this constant is quite similar to the constant 
ko of Theorem 2.2. 
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2.3 Free boundary problems 

Our proofs of Theorem 2.1 and Theorem 2.2 rely on the following two results for free 
boundary problems involving degenerate and nonlipschitz parabolic equations. We state 
them independently because they could be interesting in themselves. 

Theorem 2.5 Let x\ < x 3 < x 4 < x 2 , and let c > and T > be real numbers. 
Let g : (0, c ) — > (0, +oo) be a continuous function. Let v : (xi,x 2 ) x [0, T) — > R be a 
function such that 

(vl) v is continuous in (xi,x 2 ) x [0, T); 

(v2) v(x,t) > for every (x,t) G (x 1 ,x 2 ) x [0,T); 

(u<?) f (x, 0) > for every x G (£3,^4); 

(u^) i/ie partial derivative v x (x,t) exists for every (x,t) G (xi,x 2 ) x (0, T); 
(v5) setting 

V := {(x, t) G (x 1 ,x 2 ) x (0, T) : < v(x, t) < c } , 
we have that v G C 2,1 (V), and 

v t (x, t) = g(v(x, t))v xx (x, t) V(ar, t) G "P. 

T/jen t) > for every (x,t) G (x 3 ,x 4 ) x [0, T). 

Theorem 2.6 Let r\ < r 3 < r 4 < r 2; and /e£ Co, Ci, T 7 G 7 A be positive real numbers. 
Let g : (0, c ) — > (0, +00) fre a continuous function such that 

lim ^ = G. (2.12) 

Let f : (ri,r 2 ) x (0,T) x [— ci,ci] 2 — > R be a function such that f(r, t,0, 0) = 
uniformly in (r,t), namely 

lim sup{|/(r,f,p,g)| : (r,t,p,g) G (n,r 2 ) x (0,T) x [-a, a] 2 } = 0. (2.13) 
0-^0+ 

Lei v : (ri,r 2 ) x [0,T) — > R fre a function such that 
(vl) v is continuous in (ri,r 2 ) x [0,T); 

r(r,t) > /or even/ (r,t) G (ri,r 2 ) x [0,T); 
(v<?) f (r, 0) > /or every r G (r 3 ,r 4 ); 
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(v4) the partial derivative v r (r,t) exists for every (r,t) G (r 1 ,r 2 ) X (0, T); 
(v5) setting 

V := {{r, t) G (n, r 2 ) x (0, T) : < u(r, t) < c } , 
we have that v G C 2,1 (V), and 

v t >g(v){v rr + f(r,t,v,v r ) + A} V(r,t) eV. (2.14) 

Finally, let T> be the set defined as in (2.9) with ho := G\fA. 
Then v(r,t) > for every (r,t) G T>. 

3 Proofs 

3.1 Proof of Theorem 2.1 

We limit ourselves to prove that u x (x, t) < 1 for every (x,t) G (#3,2:4) x [0, T). The 
proof of the symmetric inequality u x (x,t) > —1 is completely analogous. 

Let us introduce some notation. Let us consider any function h G C 1 (M) which is 
nondecreasing and such that h(a) = <p'(o-) for every a G [0,1], h(a) = <f'(l) for every 
a > 1, and h(a) is constant for o < —1/2. Figure 1 shows the typical graph of such a 
function h. Note that condition (2.2) is essential for the C 1 regularity of h. 




<f'(a) h{a) g(a) 

Figure 1: Typical graph of functions <p', h, and g 

The function h, as well as the function if', is invertible as a function from (0, 1) to 
(0, We can therefore define g : (0, — > K. by setting 

g{a) := V " {h-\ V \l) - a)) Vcr G (0, (p'(l)). 

It is not difficult to see that the function g is well defined, positive, and continuous 
(but not Lipschitz continuous). In the case of the Perona-Malik equation with some 
computations one finds that g(a) = (a — a 2 ) 1 / 2 + 2(cr — a 2 ). Its graph is shown in 
Figure 1. 
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Let us finally set 



v(x, t) := </?'(!) - h(u x (x, t)) V(s, t) G (xi, x 2 ) x [0, T). 



We claim that 5 and t> satisfy the assumptions of Theorem 2.5. If we prove this 
claim, then we can conclude that v(x,t) > in (#3, #4) x [0, T). This is equivalent to 
say that h(u x (x,t)) < y?'(l), which in turn is equivalent to say that u x (x,t) < 1 in the 
same region. 

So we only need to show that v fulfils assumptions (vl) through (v5) of Theorem 2.5. 

Properties (vl) through (v3) The continuity of v easily follows from the continuity of 
h and u x . Moreover v(x, t) > in (x±, x 2 ) x [0, T) because h(a) < <f'(l) for every uel 
Due to (2.7) and the fact that <p'(cr) < <//(l) when a < 1, we have that t> satisfies (v3). 

Property (v4) It is well known that u is of class C°° where |u x (x, t)| 7^ 1 (because of 
the standard interior regularity theory for parabolic equations). Therefore the existence 
of the partial derivative v x (x,t) is trivial when \u x (x, t)\ 7^ 1. The existence of v x (x,t) 
is trivial also when u x (x, t) = —1 because h(o~) is constant for a < 1/2. 

Let us consider now a point (xo,to) with u x (xo,to) = 1. We claim that in this point 
v x (xo,to) exists and is equal to 0. Let us assume that this is not the case. Then there 
exists a sequence 5k — * such that 



Up to subsequences, we can always assume that either u x (xo + o~k,to) > 1 for every 
k G N, or u x (xq + 5k, to) < 1 for every k G N. In the first case the fraction in (3.1) is 
always 0, which is incompatible with the condition stated therein. In the second case 
the fraction in (3.1) can be rewritten as 

h(u x (x + 5 k ,t )) - h(u x (x ,t )) (p'(u x (x + 5 k ,t )) - (p'(u x (x ,t )) 



When 5k — > 0, this quotient tends to ({p'(u x )) x (x , to), and we know that this deriva- 
tive exists because u is a solution of (2.6) of class C 1 . In order to find a contradiction, 
it is enough to show that this derivative is equal to 0, and this is true because it is the 
derivative in x = x of the function x — > ip'(u x (x,t )), which attains its maximum for 
x = Xo- 

Property (v5) Let us set 



v(x + 5 k , t ) - v(x ,to)) 
5k 



> v > 



VA; G N. 



(3.1) 



V := {(x, t) G (x 1 ,x 2 ) x (0, T) : < v(x, t) < <//(!)} . 



(3.2) 
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From the properties of h it follows that < u x (x, t) < 1 in V, hence u is of class C°° 
in V. Since h and <// coincide in (0, 1), we have therefore that 

v t = -h'{u x )u xt = -h'(u x ) (<p'(u x )) xx = -<p"(u x )(h(u x )) xx = <p"(u x )(-h(u x )) xx = 

= tp"(u x )(tp'(l) - h(u x )) xx = (p"(u x )v xx 

in V . Moreover in V we can express u x in function of v as u x = h^ 1 (ip'(l) — v). In 
conclusion 

v t = (p"(u x )v xx = if" (h^ 1 (v?'(l) - v)) v xx = g(v)v xx , 
which proves that v satisfies (v5). 

3.2 Proof of Theorem 2.2 

The argument is similar to the proof of Theorem 2.1. We define h, g, and v as in that 
proof, and we claim that the assumptions of Theorem 2.6 are satisfied. 

Properties of g As in the proof of Theorem 2.1 it is quite easy to show that g : 
(0,(/9'(l)) — > (0, +oo) is a continuous function. It remains to compute the limit in 
(2.12). Since we deal with positive functions, we can square the numerator and the 
denominator. Applying the variable change r := /i _1 ((/?'(l) — a) and De L'Hopital's 
rule, we have therefore that 

, im kW]! =lim b--c-'W)-^ =lim V(r)f _ 



cr^o+ a cr^o+ a r->i- y?'(l) — ^'(t) 



which proves (2.12) with G := y/2\ip"'(l)\. 

Properties of v The proof of (vl) through (v4) is analogous to the proof of the corre- 
sponding properties in Theorem 2.1. In order to prove (v5), let us consider the set V 
defined in analogy with (3.2). As in the previous case we have that < u r (r,t) < 1 in 
this set, hence v is regular and satisfies 



v t = -h'(u r )u rt = -ip"{u r ) 



vp K)) r + — 



(3.3) 



Once again <p"(u r ) = <p" (h 1 — v)) = g(v). Moreover 

(ip'{u r )) r ip'(u r ) v r ip'(l) v 



i ii w <P'( u r) 



(p'M) rr + 
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Plugging these identities into (3.3) we obtain that 

v t = g(v) |w rr + ^ - + 1 > g{v) {v„ + f(r, t, v, v r ) + A} , 

where 

We have thus proved that i> satisfies the differential inequality (2.14) in T 3 with a 
constant A > 0, and a function f(r,t,p,q) satisfying (2.13). 

Conclusion From Theorem 2.6 we deduce that v is positive in the region V defined 
according to (2.9) with ko = Gy/A = r^ 1 \/2ip'(l)\ip"'(l)\, which in turn implies that 
u r < 1 in the same region. The proof of the symmetric inequality u r > — 1 is completely 
analogous. 

3.3 Proof of Theorem 2.3 

Let us assume that a solution exists with 

T>T , r2(r2 " ri) 



v/2yy(i)|yy"(i)|' 

Let us set 

M(t) := max{w r (r,t) : r G [r 3 ,r 5 ]} Vt G [0,T). 

Due to the first two inequalities in (2.10) we know that the subcritical region is 
nonempty at time t = 0. Applying Theorem 2.2 we have therefore that the subcritical 
region expands, and coincides with the whole interval (ri,r 2 ) as soon as 

t > = T . 

k 

In particular this means that 

M(T ) < 1. (3.4) 

On the other hand we claim that 

M(t) > M(0) - ( t~Y-t WG[0,T ]. (3.5) 
r i 

If we prove this claim, then setting t = T and exploiting the last inequality in (2.10), 
we find that 

M(T ) > M(0)-^T 



r 2 (r 2 - n ) / </?'(!) <p'(l) r 2 (r 2 -r 1 ) 



%'"(!) I r\ ^2 V '{l)\ V "\l) 
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which contradicts (3.4). 



Proof of (3.5) The argument is similar to the usual comparison principles. Setting for 
simplicity v(r,t) := u r (r,t), we have that v is a solution of 

v t = <p"(v)v rr + <p"'(v)v 2 r + ^-v r - ^ (3.6) 
in the subset of (r 1 ,r 2 ) x [0, T) where \v\ ^ 1. Let us set 

«;(*):= M(0) - e - Wg[0,T ], (3.7) 

where e > is small enough so that w(T ) > 1, hence w(t) > 1 for every t G [0, To]. We 
claim that 

M(t)>w(t) Wg[0,T ], (3.8) 

from which (3.5) follows by letting e — > + . Let us prove (3.8) by contradiction. Let us 
assume that M(t) < w(t) for some t G [0, T ], and let us set 

to := inf {t G [0,T ] : v(r,t) < w(t) Vr G [r 3 ,r 5 ]}. 

Since M(0) > w(0), we have that to > 0. Moreover, due to the continuity of v and 
w, there exists r G [r3,r 5 ] such that v (r ,to) = w(*o), an d 

f(r,i ) - w(t ) < Vre[r 3 ,r 5 }. 

Since subcritical regions don't shrink, we have that v(r 3 ,t Q ) < 1 and v(r 5 ,t ) < 1, 
while w(to) > 1. This shows in particular that tq ^ r 3 and tq ^ r$. Now we know that 
ro is a maximum point for the function r — > v(r,to) — w(to), and tq is contained in the 
open interval (r 3 ,r 5 ), hence 

f r (r ,to) — and f r r( r O)^o) < 0. (3.9) 

Let us consider now time derivatives. Since v(r ,t ) = w(to) > 1, we can use (3.6). 
Exploiting also (3.7) and (3.9) we obtain that 

v t (r , t ) - w t {t ) = <f"(v(r , t ))v rr (r , t ) - + 

> (p"(v(r ,to))v rr (r ,to). 

Since (p"(v(ro,to)) < 0, we can conclude that 

v t (r ,to) -w t (t ) > (f"(v(r ,t ))v rr (r ,t ) > 0. 

This implies that v (r , t)—w(t) > for every t in a suitable right-hand neighborhood 
of to, which contradicts the definition of to. 
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3.4 Proof of Theorem 2.4 

Let us set 

k\ := n, k 2 := 1, hi := n 3 , h 2 := —n 2 . (3.10) 

We claim that the conclusions of statement (1) and statement (2) hold true provided 
that n is large enough. 

Statement (1) Let us set for simplicity vo(x, y) := \Vuo(x, y)\ 2 . Let us assume that 

U0y(0,0)>0. (3.11) 

Then the implicit function theorem implies that the set vo(x,y) > 1 can be repre- 
sented, in a neighborhood of (0,0), as y > g(x), where g is a suitable function defined 
in a neighborhood of x — 0. Such a function satisfies 

9(01 = °- 9 (0) = 

g"(0) = Y {vQ x V yy + VQ y VQ XX - 2V Q X V QyV Q X y} , (3.12) 

where in (3.12) all partial derivatives of v are computed in (0,0). In particular g is 
convex in a neighborhood of if the right-hand side of (3.12) is positive. 
From (2.11) we have that, up to higher order terms, 

u x(x, y) = — + 2k ± x + 3/iix 2 , u 0y (x, y) = — + 2k 2 y + Sh 2 y 2 , 

hence 

v (x, y) = l + 2V2(hx + k 2 y) + (Ak\ + 3V2h^j x 2 + (4k 2 2 + 3^2) y 2 - 

All the derivatives appearing in (3.11) and (3.12) can be easily computed. It follows 
that condition (3.11) is equivalent to k 2 > 0, while g"(0) > if and only if 

8k 2 kj + 3V2 [k\h 2 + k\hx) < 0. 

Both conditions are satisfied if the values of the parameters are given by (3.10) and 
n is large enough. 
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Statement (2) Let us set for simplicity v(x,y,t) := \Vu(x, y, t)\ 2 . Thesis is proved if 
we show that 

V t = 2U X U tx + 1UyU t y > 

in the point (x, y, t) = (0, 0, 0). We can therefore deduce the value of v t (0, 0, 0) from the 
Taylor expansion of uo(x,y). 

In order to compute u t , we recall that u t = ty lx + ^/ 2y , where 



», : =*,'(k+<) 1/2 ) 



* 2 :=^(K+«;) 1/2 ) 



u,, 



With some computations we obtain that, up to higher order terms, 



(u 2 x + u 2 y ) 1/2 = 1 + V2(k lX + k 2 y) + 



fci + — fti 



x 



ko H — «2 



y 2 - 2kxk 2 xy, 



hence 



and therefore 



<f' (K + <) V2 ) = ^'(1) + (hx + k 2 y) V"(l), 



*i = ^(1) + (fez - fey) + i ^ 

<//(!) (V2kl - 3/i 2 ) + ^"'(l)fci] y 2 + (p'(l) + </'(!)) fefexy. 



1 

+ 2 



The expression for \l/ 2 is just the symmetric one. It follows that, up to higher order 
terms, 

u t (x,y,0) = (fci + fc 2 y(l) 

+ { (3^i + fci(fc 2 - 3fci)) ^'(1) + v^^"(l)A;i(A;i + A; 2 )j x 

+ { (^3/i 2 + V2k 2 (k 1 - 3k 2 )^j <p'(l) + v / 2</7 /// (l)A; 2 (A;i + A; 2 )} y, 



hence 



u t (0, 0, 0) = <//(!) {3v^(/ii + /12) + 4fciA: 2 - 6A; 2 - 6A; 2 } + 2^ /,/ (l)(fc 1 + k 2 f 



From this expression it is easy to see that v t (0, 0, 0) > if the values of the parameters 
are given by (3.10) and n is large enough. 
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3.5 Proof of Theorem 2.5 

Let x± G (#3, #4) be any point. We have to prove that 

v(x ic ,t)>0 We[0,T). (3.13) 

To this end we fix some notation. First of all we choose real numbers £5 and £ 6 such 
that £3 < £5 < x± < xq < £4. Then we consider the functions 

ip(x) := (x — x 5 )(xa — x), w(x, t) := e~ xt (5 2 il){x) + 5vp 2 (x)) , 

where A and 5 are positive parameters. We claim that when A is large enough and 8 is 
small enough we have that 

v(x,t) > w(x,t) V(x,t) G [x 5 ,x 6 ] x [0,T). (3.14) 

Since w is positive in (£ 5 ,£ 6 ) x [0, T), and £* G (xs,xq), this is enough to prove 
(3.13). In order to prove (3.14) we first establish some properties of w. 

Properties of w Let us show that w fulfils the following properties. 

(wl) w G C°°([x 5 ,x 6 } x [0,+oo)); 

(w2) w(x,t) > for every (x,t) G (x 5 ,£ 6 ) x [0, +00); 

(w3) w x (x 5 ,t) = 5 2 e~ xt ip'(x 5 ) > for every t > 0; 

(w4) w x (x 6 ,t) = 5 2 e~ xt ?p'(x 6 ) < for every t > 0; 

(w5) if 5 is small enough we have that w(x, 0) < 0) for every x G [£5, x 6 ]] 

(w6) if <5 is small enough we have that w(x, t) < c for every (x, t) G (£5, xq) x [0, +00); 

(w7) if S is small enough and A is large enough, then w satisfies 

w t (x,t) < g(w(x,t))w xx (x,t) V(x,t) G (x 5 ,x 6 ) x [0, +00). (3.15) 

Properties (wl) through (w4) easily follow from the definition of w and ip. Property 
(w5) follows from the fact that the infimum of v(x, 0) for x G [£5, £6] is strictly positive 
due to (v3). Property (w6) is almost trivial. In order to prove (w7) we recall that 
ip"(x) = —2, hence 

w t (x, t) = -Xe- Xt (5 2 i){x) + # 2 (£)) , w xx (x, t) = e- xt (25[4)'{x)} 2 - A5^{x) - 25 2 ) . 
Plugging these identities in (3.15), we are left to prove that 

-\(5ip + i) 2 ) < g(w) (2[^'] 2 -Aip- 25) V(ar, t) G (x 5 , x 6 ) x [0, +00). (3.16) 
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To this end we fix once for all two real numbers x 7 and Xg such that x 5 < x 7 < x s < 
x 6 , and 

inf {2[^\x)f -AtP(x) : x G [x 5 ,x 7 ] U [x s ,x 6 ]} > 0. 

This is possible because in the endpoints of the interval [a^a^] one nas that ip = 
and if)' 7^ 0. Now we distinguish two cases. When x G [x 5 , x 7 ] U [x 8 , x 6 ] the left-hand side 
of (3.16) is negative, while the left-hand side is positive provided that 5 is small enough, 
independently on x. When x G [a^xs] the right-hand side may be negative, but also 
the left-hand side is strictly negative because in this interval ip is bounded from below 
by a positive constant. In other words, in [a^a^] x [0, +oo) inequality (3.16) holds true 
if we choose 

f g(w(x,t)) ttibtx) + 25) , . . . . J 

We point out that the supremum is finite. This completes the proof of (3.16). 

Proof of (3.14) Let us choose positive values of 5 and A in such a way that w satisfies 
(wl) through (w7). Now we argue more or less as in the proof of the classical comparison 
results. Let us assume that (3.14) is not true, and let us set 

t := inf {t G [0, T) : 3x G [x 5 , x 6 ] such that v(x, t) — w(x, t) < 0} . 

From (w5) we have that to > 0. Moreover, from the definition of to it follows that 

v(x, t) - w(x, t) > V(x, t) G [a*, x 6 ] x [0, t ]. (3.17) 

Finally, from the continuity of v and w we deduce that there exists x G [x 5 , xq] such 
that v(xo, t ) — w(xo, t ) = 0. 

We claim that x 7^ x 5 . Indeed let us assume by contradiction that x — x 5 . Then 
w(xQ,to) = 0, hence also v(xo,t ) = 0. By (v2) it follows that x is a minimum point 
for the function x — > v(x,t ). By (v4) we have therefore that v x (xo,t ) = 0. Keeping 
(w3) into account, we deduce that 

(v - w) x (x , t ) = v x (x , t ) - w x (x , t ) < 0. (3.18) 

On the other hand, from (3.17) we know also that x = x 5 is a minimum point for 
the function x — > v(x, to)—w(x, to) restricted to the interval [x$, xq\. Since the minimum 
point is the left-hand endpoint of the interval, we deduce that (v — w) x (xo, to) > 0, which 
contradicts (3.18). 

In a symmetric way we prove that x 7^ x Q . So we are left with the case in which 
to > and x G (x 5 ,xq). In this case by (w2) and (w6) we have that < v(x ,t ) = 
w(xq, to) < Co, hence both v and w are smooth in a neighborhood of this point and fulfil 
(v5) and (w7), respectively. In particular, since xq is always a minimum point of the 
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function x — > v(x, t ) — w(x, t ), and now x is in the interior of the interval (x 5 , x 6 ), we 
have that 

v x (x ,t ) =w x (x Q ,t ) and v xx (x Q ,t Q ) > w xx (x ,t Q ). (3.19) 

Let us consider now the time derivatives. On the one hand, from (3.17) we deduce 
that (v — w) t (x ,t ) < 0. On the other hand, from (v5), (w7), and (3.19) we have that 

(v - w) t (x , t ) > g(v(x , t ))v xx (x , t ) - g(w(x , t ))w xx (x , t ) 
= g(w(x , t )) (v xx (x , t ) - w xx (x , t )) 
> 0. 

This rules out the last case and completes the proof of (3.14). 
3.6 Proof of Theorem 2.6 

The strategy is similar to the proof of Theorem 2.5. The main difference is that in this 
case we have to cope with moving domains. 

Let (r±,t+) be any point of V. We have to prove that 

u(r*,f*)>0. (3.20) 

To this end we fix some notation. First of all it is not difficult to see that there exists 
real numbers r 5 , r 6 , k such that 

| A; | < GVA, r x < r 5 + kt* < r* < r 6 + kt* < r 2 - 

Then we consider the set 

V* := {(r, t) e (n, r 2 ) x[0,T):t< U, r 5 + kt < r < r 6 + kt} C V. 

We refer to Figure 2 for a representation of the set V* (corresponding in that case 
to some k < 0) and its relation with (r*,i*) and V. Note that the slope of the slanted 
lines bounding V* is larger than the slope of the slanted lines limiting V. This is just 
because \k\ < G\[A. 

Due to this inequality, there exists e € (0, min{l, G, A/2}) such that 

\k\ < (1 - e )(G - e )VA-2eo. (3.21) 

From now on, e is a fixed positive constant. Due to (2.12) there exists also c 2 € 
(0, Co) such that 

9 {a)>{G-e )y/d WE(0,c 2 ). (3.22) 
Finally we consider the functions 

i/j{r) := (r - r 5 )(r 6 - r), w(r,t) := 5 3 ip(r - kt) + 5[tp(r - kt)f /2 , 
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r 




Figure 2: the sets V and £>* 



where 5 is a positive parameter. We claim that when 5 is small enough we have that 

v{r, t) > w{r, t) V(r, t) G V*. (3.23) 

This inequality, applied with (r,t) = (r±, i*), implies (3.20). In order to prove (3.23) 
we first establish some properties of to. 

Properties ofw Let Int(X>*) denote the set of points (r, t) G Z>* with r 5 +kt < r < r 6 +kt 
Let us show that w fulfils the following properties. 

(wi) tueCH^nC 00 ^)); 

(w2) w(r,t) > for every (r, t) G Int(Pj; 

(w3) m r (r 5 + kt,t) = 5 3 4}'(r 5 ) > for every t G [0,4]; 

(w4) w r (r 6 + kt,t) = 5 3 -0'(r 6 ) < for every t G [0,t*]; 

(w5) if 5 is small enough we have that w(r, 0) < v(r, 0) for every r G [r 5 , r 6 ]; 

(w6) if S is small enough we have that w(r, t) < c 2 for every (r, £) G P*; 

(w7) if 5 is small enough, then w satisfies 

w t < g(w) {w„. + f(r, t, w, w r ) + A} V(r, t) G Int(X\). (3.24) 

Properties (wl) through (w4) easily follow from the definition of w and ip. Property 
(w5) follows from the fact that the infimum of v(r, 0) for r G [r 5 , r 6 ] is strictly positive 
due to (v3). Property (w6) is almost trivial. In order to prove (w7) we recall that 
ip"(r) = —2, hence (for simplicity we set y := r — kt, and we observe that y G [r$, r^]) 

w t (r,t) = -k5^'(y) - hd [^(y)} 1/2 f(y), 
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w r (r,t) = ^(y) + ^(y)] 1/2 ^(y), 

w rr (r, t) = -25 3 - 35 [4>{y)] l/2 + \d h%)]~ 1/2 W{y)f . 

When 5 — > + we have that w and w r tend to zero uniformly in T>*. Thanks to (2.13) 
we have therefore that 

\f(r,t,w(r,t),w r (r,t))\ < e V(r,t) G D*, 

provided that S is small enough. In an analogous way we have also that 

w rr (r,t)> 3 -5myT 1/2 W(y)] 2 -eo 
provided that S is small enough. It follows that 

w rr + f(r, t,w,w r ) + A> A - 2e + [ip{y)]~ 1/2 W{y)f > 
in Int (£>*). Moreover, from (3.22), (w2) and (w6) we have that 

g(w(r,t)) > (G-e )y/w(r,t) V(r,f) G Int(P*), 

and in conclusion 

g(w) {w rr + f(r, t, w, w r ) + A} > (G - e )Vw I A-2e + -5 ^- 

Therefore inequality (3.24) is proved if we show that 

8*\ktf\ + p\kip'\^p< (G - E )^5^ + 5[^/ 2 {A - 2s + > ( 3 - 25 ) 

where the argument of ip and ip' is any y G (r 5 ,r 6 ). Let us consider the right-hand side 
of (3.25) multiplied by (1 — e ). Applying the inequality between the arithmetic mean 
and the geometric mean, and recalling (3.21), we obtain that 

(l-e )- (right-hand side) > (1 - e )(G - e )^5[^/ 2 • |a - 2e + 

> (l-eo)(G-eo) V / ^^-2[(^-26o)-^^| 
= (1 - e )(G - e )y/A - 2e ■ 8\^'\ ^ ■ Vz 



1/2 



2 
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In order to prove (3.25) it is therefore enough to show that 

5 3 N'| < e (G - e )^5^ + 5[^/^A - 2e + , 
which in turn is true if we show that 

5 3/ Wl < e (G - e ) {^4 - 2e + §<^} ■ 
To this end, we fix once for all two real numbers r-; and r 8 such that 

r 5 <r 7 < — - — < r 8 < r 6 . 

When y G (r 5 , r 7 ] U [r 8 , r 6 ) we have that ^'(y)! is bounded from below by a positive 
constant. Therefore 

5*' 2 \W\y)\<e Q {G-e Q )- A 5W{y)¥ 

provided that 5 is small enough. When y G [r7,r 8 ] we have that ^(y) is bounded from 
below by a positive constant, hence in this case 

5 3/2 \k4>'(y)\ < e (G - e )(A - 2e Q )^y) 

provided that 5 is small enough. This completes the proof of (3.25) and shows that w 
satisfies (w7) whenever 5 is small enough. 



Proof of (3.23) The argument is analogous to the proof of the corresponding inequality 
in Theorem 2.5. Let us choose a positive value of 5 in such a way that w satisfies (wl) 
through (w7). Let us assume that (3.23) is not true, and let us set 

t := inf {t G [0, t*\ : 3r G [r 5 + kt, r 6 + kt] such that v(r, t) - w(r, t) < 0} . 

From (w5) we have that to > 0. Moreover, from the definition of to we have that 

v(r,t)-w(r,t)>0 V(r,t) G £>* with t < t . (3.26) 

Finally, due to the continuity of v and w, we deduce also that there exists r G 
[r 5 + fct , r 6 + kt ] such that v(r , t ) - w(r , t ) = 0. 

We claim that r ^ r 5 + kt . Indeed let us assume by contradiction that r = r 5 + /cto- 
Then w(ro,to) = 0, hence also i>(ro,to) = 0. By (v2) it follows that r$ is a minimum 
point for the function r — > t>(r, t ). By (v4) we have therefore that f r (r ,t ) = 0. 
Keeping (w3) into account, we deduce that 

(v - w) r (r , t ) = v r (r , t ) - w r (r , t ) < 0. (3.27) 
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On the other hand, from (3.26) we know also that r = r 5 + kt is a minimum 
point for the function r — > t> (r, t ) — w(r, t ) restricted to the interval [r 5 + /ct ? r 6 + ^o]- 
Since the minimum point is the left-hand endpoint of the interval, we deduce that 
(v — w) r (ro,to) > 0, which contradicts (3.27). 

In a symmetric way we prove that r$ ^ r§ + Hq. So we are left with the case in 
which t > and r G (r 5 + kt ,r 6 + kt ). In this case < t>(r ,to) — w ( r o?^o) < c 2, 
hence both v and w are smooth in a neighborhood of this point and fulfil (v5) and 
(w7), respectively. In particular, since r is always a minimum point of the function 
r — > f(r, t ) — w^^o), and now r is in the interior of the interval (r 5 + kto,r 6 + /ct ), 
we have that 

v r (r ,t ) = w r (r ,t ) and v rr (r ,t ) > w rr (r ,t ). (3.28) 

Let us consider now the time derivatives. On the one hand, from (3.26) we deduce 
that (v — w) t (r ,t ) < 0. On the other hand, from (v5), (w7), and (3.28), in the point 
(r ,to) we have that 

(v-w)t > g(v) {v rr + /(r, t, v, v r ) + A} - g(w) {w rr + f(r, t, w, w r ) + A} 
= g(w) {v rr + f\r, t, w, w r ) + A} - g(w) {w rr + /(r, t, w, w r ) + A} 
= g{w) (v rr — w rr ) > 0. 

This rules out the last case and completes the proof of (3.23). 
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